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C D ' Well above the critical temperature hot QCD is described by 3d electrostatic QCD 

C D ' with gauge coupling gE and Debye mass ttie ■ We integrate out the Debye scales to 

£\J i two loop accuracy and find for the gauge coupling in the resulting magnetostatic 
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J"""-. ' 1 Introduction 

Notable progress 1 in the standard model at high T is due to the systematic separation of 
(f) ' perturbative scales like the temperature, the Debye mass gT and the non-perturbative scale 

g 2 T. It is the latter that is three dimensional and can be treated numerically on the lattice 
and has given us a wealth of information on the plasma state of the standard model and QCD 
itself. In this note we will be concerned with QCD, but we will admit for N instead of three 
colours. 

For small gauge coupling g one can integrate out the integer Kaluza-Klein modes 2nnT 
and obtain a static effective action, Se, with a running coupling gE(T). This action is three 
dimensional, and its degrees of freedom are three dimensional Yang-Mills and a massive adjoint 
Higgs with as mass the Debye scale gT. The scale g 2 T appears as the coupling in the three 
Qh! dimensional Yang- Mills action. The Debye scale can be integrated out when g « 1. We are 

then left with the magnetic action Sm- It is the three dimensional Yang-Mills theory which 
, describes physics at the magnetic scales g 2 T '. Its coupling g\ L is a function of the parameters 

in Se and can be calculated perturbatively. This has up to now been done to one loop order 4 . 
In this letter we report on the computation of the two loop effects. 



2 Motivation 

Our motivation stems from the need for accuracy. More precisely, integrating out effects of the 
K-K modes one obtains from the original action of QCD the superrenormalizable action Se- 

C E = Tr(D(A)Ao) 2 +m 2 E TrA 2 + \E(Tr(A 2 Q )) 2 + 

+ \E(Tr{A,f~ l -{TrAl) 2 )+ l -TrF 2 +5C E . (1) 

This action density describes the physics of the QCD plasma down to temperatures T ~ 2T C , 
including non-perturbative effects from the magnetic density TrF 2 ^ . These have been calculated 
by lattice methods 2 . The terms neglected, 6Ce, introduce an error of 0(g 4 ) 3 . So the parameters 
in this action have all evaluated up to this order. 

The magnetic action takes the form 



C M = -TrFZ + 6£ M (2) 



1 
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with a magnetostatic gauge coupling g 2 M . 

Now the neglected terms introduce an error 3 of 0(g 3 ), and hence the magnetic coupling 
has to be computed to 0(g 2 ) accuracy. The calculation is reported on in the next section. 
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3 Renormalization of the magnetic gauge coupling 



The basic idea behind the effective actions eqns (1) and (2) is that one can compute with both 
in the region of momenta p ~ g 2 T . To know what the parameters of the latter are in terms of 
those of the former requires computing two-point functions, three point funtions etc. in both 
theories and match them. In the matching the diagrams of the pure 3d Yang-Mills theory drop 
out. 

Here we will follow a well-known shortcut 5 by introducing a background field Bi in Ce- 

A = —B + Q 

9E 

A = 9eQo- (3) 
We calculate the fluctuations around the background in a path integral: 

exp -A-S M {B) = [ DQ DQiexp (-S E - ^Tr(AQi) 2 )- (4) 

We added a general background gauge term. The resulting action Sm(B) is gauge invariant 
to all loop orders and the renormalization of the coupling is identified from the background 
field two point function at a momentum p = 0(g 2 T). This momentum is the infrared cut-off 
in computing the r.h.s. of eq.(4). With dimensional regularization one finds in d dimensions, 
dropping the pure Yang-Mills diagrams as mentioned before: 

exp (--Ls M (B)) = exp{-\s M {B))(l + {F{ r + F? + ...)S M (B)). (5) 
9m 9e 

Here the F\ r are the transverse parts of the background two-point functions Fi as shown for 
the two loop case in fig.(l). 

Let the sum of all Feynman diagrams for the two point function of the background field 
with i loops be Fj. Then we can write: 

Fi = F\ r {5 lmP 2 - p lPm ) + Ft PlPm (6) 

The longitudinal part is zero because of gauge invariance. It is borne out by explicit 
calculation. Transversality is true for all d and values of the parameters. 

The Ff r are still depending on the parameters e = gs, rriE, £, the regularization scale 
p and the momentum p. We are interested in the limits d = 3 and p = where p = In 
that limit we expect because of the superrenormalizabilty the UV poles and the j2 dependence 
to cancel. Also the £ dependence should disappear. And so should the IR effects in the guise 
of inverse powers of the momentum p. And indeed they do by explicit calculation, as the table 
in the next section shows. 

This leaves us with the relation, using eq.(5): 

4- = 4" - F t - ( 7 ) 

9m 9e 

with 



This is the main result. 



2 F tr ... 19 ( 9%N 2 



2 




massive propagator 
massless propagator 
background field 
one loop renormalisation 

Figure 1: two-loop Fcynman graphs for the background two-point function 



4 Details of the calculation 



The two loop diagrams involving at least one massive propagator are shown in fig.(l). Also 
shown (graphs 15 and 16) are the insertions, discussed in ref.f ). They are vital for the gauge 
parameter independence of our result. FORM 7 was used for algebraic manipulations and 
scalarisation of integrals with reducible numerators. The program TARCER 6 served to scalarize 
those integrals with irreducible denominators. At that point the result F2 is expressed in terms 
of eight scalar integrals. For d = 3 , the values of these scalar integrals are computed in ref.f). 
For the expansion of these integrals up to e, methods as in ref.f) were used and results were 
checked with ref.f). 

First we checked the transversality, i.e. F% = 0. The reader can find the result for F| r in 
table below. Individual graphs have UV and IR divergencencies that do cancel when summed. 
The finite part is indeed gauge parameter independent, though the physically irrelevant 0(p 2 ) 
terms are not. 



Graph no 



v 2 n r 



1 

2 
3 
4 
5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 



sum 



p 2 (29+24Q , l+g+ £ (2+£-4(l+C) log(2))+4e(l+g) log(£) 

92167T 2 128£7T 2 

p 2 (20+6Q _ 3(2+g-2e(-l+g(-l+log(4))+log(16))+4e(2+g)log(£))) 

f (-43+12Q 4+ £ (l+g-161og(2)) + 16elog(^) 
4608tt 2 1 " 64^ 

2-_i _ p(15+!) ^ p 2 (13+2g) _ 3(l+g) 2 +46(-8+Oeiog(£ r )+4e(3+2g(7+g))log(^) 



32p~7T 



1536tt 



+ 3072tt 2 
p(15+C(4+Q) 



p 2 (23+12g) 



1536tt 



-3 
(64e7r a ) 



4608tH 1 
. 22-24(l+log(^)) 



768e7r 2 £ 
2e+g +4e? log(^) 

128CT 2 



1287T 2 

6(2+|)+e(12(l+g) -p 2 (2+g))+24e(2+$)log(; 

512£7T 2 





l+2g 
256e7r 2 £ 



+ 



-2g-4(l+24)-4(-l-2g)(l+log(^| r )) 
256tt 2 £ 



16+8g 
256p7r 








p(l+g)(4+g) 
1536tt 
P(2+C) 
768tt 



-19p 2 
4608tt 2 



5 Conclusions 

Our main result, eq.(9), shows that the smallness of the corrections to the magnetic coupling 
does persist in two loop order. In fact, at 2T C the coupling for 3 colours equals 2 g% = 2.7 and 
the 2 loop correction is about a third of the one loop correction (itself about 3 percent). 

Our result is of importance in analyzing the purely magnetic quantities, as the spatial 
Wilson loop, and the magnetic mass. In particular it is crucial in connecting the lattice results 
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from the magnetic action to those obtained from the electric action, and ultimately to those 
of four dimensional simulations. This will be done in a future publication 10 . I thank Chris 
Korthals Altes for his help throughout my work. I acknowledge the help of Mikko Laine and 
York Schroder for very useful and stimulating advice. I thank the MENESR for financial 
support. 
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